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The processes of ionization and energy transfer in a quantum system composed of two distant H atoms with an initial internuclear 
separation of 100 atomic units (5.29 nm) have been studied by the numerical solution of the time-dependent Schrodinger equation 
,beyond the Bom-Oppenheimer approximation. Thereby it has been assumed that only one of the two H atoms was excited by 
■temporally and spatially shaped laser pulses at various laser carrier frequencies. The quantum dynamics of the extended H-H 
system, which was taken to be initially either in an unentangled or an entangled ground state, has been explored within a linear 
three-dimensional model, including two z coordinates of the electrons and the internuclear distance R. An efficient energy transfer 
|from the laser-excited H atom (atom A) to the other H atom (atom B) and the ionization of the latter have been found. It has been 
shown that the physical mechanisms of the energy transfer as well as of the ionization of atom B are the Coulomb attraction of 
the laser driven electron of atom A by the proton of atom B and a short-range Coulomb repulsion of the two electrons when their 
wave functions strongly overlap in the domain of atom B. 
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Introduction 



Long-range energy transfer among distant quantum systems 



> 



^ ,has attracted considerable interest recently in theory and ex- 
Q^periment.— ~— Extended quantum systems such as those com- 
posed of two atoms separated by several nanometers are also 
referred to as 'extreme quantum systems'. ^ The energy trans- 
|fer from an excited atom to its neighbor, the so-called in- 
teratomic Coulombic decay (ICD), has been first studied by 
,Cederbaum et al. for molecular clusters.— ICD is well es- 
■tablished experimentally for inner- valence excitation of many 
electron systems. In recent work^, ICD was demonstrated 
experimentally for a helium dimer at interatomic distances up 
|to about 12 a.u. Since helium atoms have no inner- valence 
electrons, a different type of ICD has been suggested for this 
• • case which does not require an overlap of the electronic wave 
^ functions.— 

The present work is addressed to an extended quantum sys- 
tem composed of two H atoms, to be referred to as atoms A 
.'and B. The atoms are assumed to be initially separated by an 
internuclear distance of R — 100 a.u. which, although chosen 
somewhat arbitrary, will highlight the particular properties of 
the laser driven dynamics of such systems. The extended H-H 
system interacts with laser pulses having narrow spatial en- 
velopes, such that practically only atom A is excited by the 
laser field directly, while atom B is almost not affected. The 
laser field is linearly polarized along the z-axis, and the pro- 
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tons and electrons are assumed to move along the polariza- 
tion direction of the laser field only. Three different initial 
conditions will be considered, which comprises one being of 
atomic origin where the (unentangled) electronic wave func- 
tion is given by the direct product of the respective atomic 
wave functions, and two corresponding to a molecular origin 
where the (entangled) electronic wave function is either sym- 
metrized (singlet) or antisymmetrized (triplet). 

In our recent work-, energy transfer and ionization in this 
extended H-H system excited by spatially shaped laser pulses 
(both narrow and Gaussian type shape) have been studied at 
the laser carrier frequency of CO = LO a.u., corresponding to 
the ground-state energy of H-H at R = 100 a.u.. Both unen- 
tangled atomic states and entangled molecular singlet states 
were used as the initial states in ref. 9. In the present work, 
we extend this study to a wide range of laser carrier frequen- 
cies, including co = 0.5 a.u. corresponding to the ground-state 
energy of an individual H atom. Further, we will include the 
case of an initial triplet state. As far the spatial shape of the 
laser field is concerned we will consider the case of a narrow 
shape only, which excites dominantly one atom. Therefore, 
such a pulse is most suitable for investigating energy transfer 
and ionization in this extended quantum systems. For a dis- 
cussion of the effect of a broad spatial envelope we refer to 
ref. 9. 

The relative simplicity of the H-H model used here, i.e. 
one-dimensional electrons (coordinates zi and zi) and protons 
(coordinate R), makes it possible to treat the long-range elec- 
tronic motion explicitly together with the nuclear motion and 
to reveal the role played by the overlap of the electronic wave 
functions and by the electron-electron Coulombic repulsion. 
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However, effects similar to those presented below can be an- 
ticipated to occur in more complex extended quantum systems 
as well. Such systems can be in the bound state, as a he- 
lium dimer^ for example, or represent spatially well separated 
atoms as those studied in the present work. 

The paper is organized as follows. The model, equations 
of motion and techniques are described in section |2] Results 
obtained for energy transfer and ionization of H-H in a wide 
range of carrier frequencies of the applied temporally and spa- 
tially shaped laser pulses and the laser-driven dynamics of H- 
H at several carrier frequencies are presented in section |3] The 
results obtained are summarized in the concluding section |4] 

2 Theoretical Model and Computational Meth- 
ods 

2.1 Time-dependent Schrodinger Equation 

Within the linear three-dimensional model used in this paper 
to describe the H-H system excited by temporally and spatially 
shaped laser fields, the total Hamiltonian Hi consists of two 
parts, 



(1) 



where Hs{R,zuZ2) represents the H-H system and 
HsF{zi,Z2,t) describes its interaction with the laser field. 
The applied laser field is linearly polarized along the z-axis, 
the nuclear and the electronic motion are restricted to the 
polarization direction of the laser electric field. Accordingly, 
the two z coordinates of electrons, zi and Z2, measured with 
respect to the nuclear center of mass, are treated explicitly 
together with the internuclear distance R. A similar model has 
been used previously, e.g. in ref. (lOl for the H2 molecule. 

The total non-Bom-Oppenheimer Hamiltonian of the H-H 
system (employing atomic units) reads 



HsiR,zuZ2) 



m. 
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where nip is the proton mass, /jg = 2mp/(2mp + 1) is the elec- 
tron reduced mass, and non-diagonal mass-polarization terms 
are neglected. The Coulomb potentials in eqn ^ are 



where k = 1,2, and the regularization parameters, a = 0.1 x 
10^^ and P = 1.995, have been chosen such as to repro- 
duce the ground-state (GS) energy of the H-H system at R — 
100 a.u. (Egs^ = -1.0 a.u.). 

The interaction of the H-H system with the laser field, 
HsF{zi,Z2,t), is treated semiclassically within the electric 
dipole approximation via the vector potential. 



' . X 1 dA(t) , 2 

c dt ^fi 



where Y = (1 +2mp)^^ A{t) is the vector potential, c is the 
speed of light, and is the spatial-shape function, or spatial 
envelope, of the laser field. The following, sin^-type, narrow 
spatial envelope of the laser pulse acting practically only on 
atom A was used for z = z\ and z = Z2'- 



F{z) 



sm 



Ti{z-Za) 
Zh - Za 



Za<Z< Zb, 



(5) 



where Za — —60 a.u., Zb — —40 a.u., and F{z) ~ otherwise. 
The vector potential A (f) in eqn (HI is chosen as follows: 



A(f) = — 'Eosin (m/fp)cos((af + (|)), 



(6) 



where To is the amplitude, tp is the pulse duration, co is the 
laser carrier frequency, and (|) is the carrier-envelope phase. 
The definition of the system-field interaction via the vector po- 
tential by eqn (|4|i, suggested in ref. IJJJ, ensures that the electric 
field has a vanishing direct-current component and satisfies 
Maxwell's equations in the propagation region. 

It is suitable to define, on the basis of eqns (|4]i and the 
local effective-field amplitude for atom A as follows: 



■E^ = -EoFizA), 



(7) 



where za = —50 a.u. (atom B is located initially at zb = 
50 a.u.). The respective time-dependent electric field acting 
on atom A reads 



K 



■ sin(2m/fp)cos(cof 



(8) 



The first term in eqn (|8]l corresponds to a laser pulse with a 
sin^-type temporal envelope, and the second term appears due 
to the finite pulse duration. ^'''^ 

The 3D time-dependent Schrodinger equation. 



^PpW = 4. ^ee(zi,Z2) 



1 



R 



\/{zi ~Z2)^ + a 



VepiZk,R) = 



v/(zfc-^/2)2 + p ^(z,+R/2)^ + ^ 



, (3) 



i-'f^ [Hs{R,ZUZ2)+HsF{zuZ2,t)]'i', 



(9) 



has been solved numerically with the propagation technique 
adapted from refs. [l3| and LL4j for both electron and proton mo- 
tion. In particular, calculations for the electron motion have 
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been performed by using 200-point non-equidistant grids for 
the Hermite polynomials and corresponding weights for the 
numerical integration on the interval (— cxi^oo) for the zi and Z2 
coordinates. For the nuclear coordinate/?, a 256-point equidis- 
tant grid has been used on the interval [75 a.u., 125 a.u.]. The 
time-step of the propagation was At = 0.021 a.u.. 

The wave function of the initial state has been obtained by 
the numerical propagation of the equation of motion, eqn (|9), 
in the imaginary time, with the laser field being switched off 
(£o — 0). Three cases will be considered, which correspond 
to different preparation conditions and are either of atomic or 
molecular origin (see next subsections). Note before proceed- 
ing that the case of unentangled atomic states is more likely 
to be realized experimentally. For example, at a gas pressure 
of 1 atm., the interatomic distance is about 100 a.u.. In con- 
trast entangled molecular states require additional steps such 
as dissociation of H2 molecules or generation of entanglement 
between two individual H atoms. 

The wave function analysis will be performed on the basis 
of the electron probabilities P{zi) and P{z2) defined as fol- 
lows: 

(10) 



P(zl/2)= / dR / dZ2/l\'i'{R,Zl,Z2) 



for electron ei/2 initially belonging to atom A/B with proton 
Pa/b- These electron probabilities give the overall probabil- 
ity to find an electron at a specified point of the z-axis at any 
position of the other electron and at any internuclear distance. 

2.2 Product Initial State 

In this case, which is termed to be of atomic origin (two distant 
H atoms), the spatial part of the initial unentangled ground- 
state wave function of H-H used in the imaginary time prop- 
agation was defined by the non-symmetrized Heitler-London 
wave function as follows: 

1'(/;,Zl,Z2,f =0)=^15a(zi)^1Sb(Z2)^gW, (11) 



where 'Pis^(zi) = e 



at ZA = -50 a.u., 'Pissfe) = 



g-k2-JBl at zb = 50 a.u., and ^g(^) is a Gaussian function 
of unit width and centered at R= 100 a.u.. 

Upon excitation of the extended H-H system by the laser 
field, the electronic wave functions of its atomic A and B parts 
may overlap. In order to study the energy transfer from atom A 
to atom B, the respective 'atomic' energies, £^{1} and £^{1), 
have been defined on the basis of eqns. © and ^ as follows: 



£a/b(0-(^(0 



2OTp dR^ 



\vp,{R)+^-V,,{z,,Z2) 



It is seen that the kinetic energy of nuclear motion and the po- 
tential energies of the proton-proton and the electron-electron 
interaction are assumed to be equally shared between atoms A 
and B. The kinetic energy of electron ei and the energy of 
Coulombic interaction of both electrons ei and 62 with proton 
PA are entirely assigned to atom A. Similarly, the kinetic en- 
ergy of electron 62 and the energy of Coulombic interaction of 
both electrons with proton pQ are entirely assigned to atom B. 
The sum of 'atomic' energies gives the correct total energy of 
the entire H-H system. 

The ionization probabilities for atoms A and B have 
been calculated from the time- and space-integrated outgoing 
fluxes, separately, for the positive and the negative directions 
of the zi and Z2 axes at zi 2 — ±91 a.u.. We calculated four ion- 
ization probabilities: Ia(zi = —91 a.u.) and Ia(z2 — —91 a.u.) 
for atom A, and Ib(zi =91 a.u.) and Ib(z2 =91 a.u.) for atom 
B. At the outer ends of the z-grids, absorbing boundaries have 
been provided by imaginary smooth optical potentials adapted 
from that designed in ref. 15. Similar optical potentials have 
also been included for the R-axis, but in practice the wave- 
packet never approached the outer ends of the R-gnd. 

2.3 Entangled Initial States 

The initial states of molecular H-H systems, singlet and triplet, 
are entangled by spin exchange and represented by the Heitler- 
London symmetrized and antisymmetrized products of atomic 
wave functions, respectively. The spatial parts of the initial en- 
tangled wave functions for singlet (S) and triplet (T) molecular 
electronic states are given by 

^{R,Zl,Z2,t = 0)s/T = [^I^aCzO'PiSb iZ2) 



±^lS^izi)^lS^{Z2)]^G{R), 



(13) 



where ^isab(zi,2) are defined by eqn (10) and ^g(^) is a 
proton Gaussian function centered at R = 100 a.u.. Imagi- 
nary time propagations have been performed with the initial 
wave functions ( fT3] l and the system Hamiltonian of eqn (2). 
The energies of the singlet and triplet molecular states of the 
extended H-H system at the large internuclear separation of 
R = 100 a.u. are equal to each other and to the energy of the 
unentangled atomic initial state. 

For the entangled initial states (T3[ . the 'atomic' energies 
£'A(f) and Es{t) are defined on the basis of eqns. (2) and (3) 
as follows: 
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It is seen that the kinetic energies of nuclear and electronic 
motion as well as the potential energies of the proton-proton 
and the electron-electron Coulombic interaction are assumed 
to be equally shared between atoms A and B. The energy of 
Coulombic interaction of both electrons with proton pp^ is as- 
signed to atom A, while the energy of Coulombic interaction 
of both electrons with proton is assigned to atom B. The 
sum of these 'atomic' energies always gives the correct total 
energy of the entire H-H system. The choice of the electron ki- 
netic energies in eqn (fT4l i is different from that used in eqn (12) 
and corresponds to the initial electron probabilities P{zi) and 
P{z2) in the entangled molecular states. Indeed, both electrons 
e\ and 62 are localized with the 50% probability in the vicinity 
of proton Pa of atom A and when the extended H-H system is 
excited by the narrowly shaped laser pulse of Fig. 1(b), both 
electrons give rise to the 'atomic' energy Ep^{t) with the same 
probability. 

Our numerical simulations have shown that the ionization 
probabilities of electrons e\ and e\ excited from the entan- 
gled initial states, singlet and triplet, are almost identical to 
each other both for the positive and the negative directions 
of the z-axes: Ia(zi — — 91a.u.) — Ia(z2 = — 91a.u.) and 
Ib(zi = 91 a. u.) = Ib(z2 = 91 a. u.). This is the consequence 
of the symmetry of the entangled initial wave functions (cf. 
ref. It is therefore reasonable to define the total ionization 
probabilities for atoms A and B as follows: 



Total 



TTotal 



= Ia(zi = -91 a.u. 
Ib(zi =91a.u. 
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Fig. 1 The unentangled initial state of the extended H-H system and 
a narrow spatial shape of the applied laser pulse, (a) - electron prob- 
abilities P{z\) and P{z2) defined by eqn Go); (b) - a narrow spatial 
envelope of the laser pulse, defined by eqn (|5}, which dominantly 
excites only electron e\ in the domain of atom A. 



3 Numerical Results 
3.1 Product Initial State 

The probability distribution for the initial state of the extended 
H-H system representing two distant H atoms, A and B, is 
shown in Fig. [TJa)- The narrow spatial envelope of the ap- 
plied laser pulses is presented in Fig. [Tib) to illustrate the local 
effective-field strengths acting on atoms A and B. We clearly 
see from Fig. [T]that the narrowly shaped laser pulse excites 
practically only electron e\ (coordinate zi) in the domain of 
atom A, because the probability to find electron 62 (coordinate 
zi) at Z2 = —50 a.u. is extremely small. 

The electric field of the laser pulses is defined by eqns (|7]i 
and (HJ, where we choose the amplitude 'Eq = 0.02 a.u. (cor- 
responding to an intensity of /o = 1.4 x lO'^ W/cm^) and the 
pulse duration tp = 5 fs. The laser carrier frequency co is 
varied in a wide range covering, in particular, (O — 0.5 a.u. 
(corresponding to the ground-state energy of an H atom) and 
CO = 1 .0 a.u. (corresponding to the ground-state energy of the 
extended H-H system). In passing we note that with the laser 



pulse duration being fp = 5 fs the number of optical cycles dur- 
ing the pulse, A^^ = a)fp/27t, is about 33 at CO = 1.0 a.u., and 
therefore, the carrier envelope phase (CEP) of the laser field, 
(|), is not important. In contrast, at CO < 0.5 a.u. the number 
of optical cycles is A^^ < 15, and CEP may play an impor- 
tant role. However, postponing a detailed study of the role 
of CEP-effects, in the present work we set the carrier phase 
equal to zero [(|) = in eqns. (|6]l and ([8])]. 

The 'atomic' energies and E-r, eqn (fT2l l. and 'atomic' 
ionizations Ia(zi = —91 a.u.) and Ib(z2 — 91 a.u.) as ob- 
tained from the numerical solution of the time-dependent 
Schrodinger equation for various laser frequencies are shown 
in Fig. |2] The 'atomic' energies and ionizations, presented in 
Fig. m are calculated at f = 10 fs, which corresponds to the 
excitation of H-H by a 5 fs laser pulse and a 5 fs free evolution 
of excited H-H after the end of the laser pulse. Notice that 
in general a small/large 'atomic' energy, E-q, and ionization, 
Ib, of atom B for an excitation of atom A by the laser pulse 
with a narrow spatial envelope [Fig. [Hb)] would evidence a 
weak/strong energy transfer from A to B. 



4 I Submitted to Phys. Chem. Chem. Phvs..1-FiTI 



o 




0.2 0.4 0.6 0.8 1 1.2 1.4 

FREQUENCY (a.u.) 

(b) 



Ia(z,=-91 a.u.) 




04 0.6 0.8 
FREQUENCY (a.u.) 



Fig. 2 Excitation of tlie extended H-H system by laser pulses with the 
narrow spatial envelope of eqn (|5) at various laser carrier frequencies: 
'atomic' energies and 'atomic' ionizations at f = 10 fs (5 fs excitation 
by the laser pulse followed by 5 fs free evolution), (a) - 'atomic' 
energies £a and defined by eqn l ll2t : (b) - 'atomic' ionizations in 
positive and negative directions of the z-axis. Note that the ionization 
probabilities Ia(zi = 91a. u.) and Ib(z2 = — 91a.u.) are of the order 
of lO^'* and not included here. 



The following observations can be made from Figs.|2la) and 
(b): 

(i) The frequency dependence of both 'atomic' energies 
[Fig- 12a)] and, especially, 'atomic' ionizations [Fig. |2b)] 
can be divided into two domains: a low-frequency domain 
(CO < 0.6 a.u.) with pronounced changes and a high-frequency 
domain (co > 0.6 a.u.) with smooth changes. 

(ii) In the low-frequency domain, including the resonant fre- 
quency CO = 0.5 a.u., which corresponds to the one-photon res- 
onant excitation of an isolated H atom from its ground state, 
the 'atomic' energies Ea and [see Fig.|2a)] first increase 
with CO to reach their maxima at CO = 0.3 a.u.. At larger fre- 
quencies, CO > 0.3 a.u., Ea decreases more rapid than Eq such 
that their numerical values approach each other towards the 
end of this interval. Overall, however, the 'atomic' energy E^ 
is substantially smaller than Ea in the low-frequency domain. 

(iii) In the low-frequency domain (co < 0.6 a.u.), 'atomic' 
ionizations, Ia(zi = —91 a.u.) for atom A and Ib(z2 = 91 a.u.) 



03=0.4 a.u. ,../ ^"i 




5 10 
TIME (fs) 



Fig. 3 Time-dependent deviation Az2{t) of the coordinate expecta- 
tion value of electron £2 of atom B [Eq. ( I16H at three laser carrier 
frequencies, co =0.4, 0.5 and 0.6 a.u., of the 5 fs laser pulse with the 
narrow spatial envelope of Eq. (|5) which excites only electron e\ of 
atom A of the extended H-H system. 



for atom B, demonstrate a very sharp and irregular dependence 
on the laser carrier frequency CO [Fig.|2jb)]. Indeed, at CO = 
0.1 and 0.2 a.u., the ionization probability of atom B is large 
compared with that of atom A, implying an efficient transfer 
of electronic energy from atom A to atom B. In contrast, at 
CO = 0.3 a.u., the ionization probability Ia is more than four 
times larger than Ib. At CO = 0.4 a.u. the ionization of atom B 
dominates again, while at CO = 0.5 a.u. ionization of atom A is 
almost nine times larger than ionization of atom B, implying 
that energy transfer from atom A to atom B is not efficient at 
a resonant excitation of atom A from its ground state. 

(iv) Finally, in the high-frequency domain (co > 0.6 a.u.), 
both 'atomic' energies and 'atomic' ionizations decrease 
smoothly when CO increases, with Ea being close to Eb, and 
Ia being close to Ib. Below we shall pay attention to the low- 
frequency domain (co < 0.6 a.u.) only, while the behavior in 
the high-frequency domain will be addressed elsewhere. 

The most striking observation in the low-frequency domain 
is that while the 'atomic' energies behave quite smoothly the 
ionization probabilities show rapid changes in particular when 
sweeping through the isolated atom resonance at CO =0.5 a.u. 
This behavior is reminiscent of the oscillatory dependence 
of the ionization probability reported for single H-atoms, al- 
though for somewhat stronger and shorter laser fields in ref. 
171 In the present case the reason of the behavior of ionization 
probability, e.g. Ib, on the laser carrier frequency CO can be 
understood from the time-dependent deviation of the coordi- 
nate expectation value for electron 62 of atom B from its initial 
position at f = 0, i.e. 



(16) 



The time-dependent deviation Az2{t) is presented in Fig. [3] 
on the timescale of 15 fs for three laser carrier frequencies: 
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Fig. 4 Excitation of the extended H-H system by tlie laser pulse with the narrow spatial envelope of eqn (|5) acting only on atom A in the 
vicinity of zh = —50 a.u.: an intuitive case (left panel, to =0.5 a.u.) and a counter-intuitive case (right panel, to =0.4 a.u.). (a) and (d) - 'atomic' 
energies £a and £b; (b) and (e) - 'atomic' ionizations, Ia(zi = —91 a.u.) for atom A and Ib(z2 = 91 a.u.) for atom B [ionization probabilities 
Ia(zi = 91 a.u.) and Ib(z2 = —91 a.u.) are of the order of 10^^ and are not shown]; (c) and (f) - laser fields acting on atom A. 



CO =0.4, 0.5 and 0.6 a.u.. From Fig. [3] we see that the dy- 
namics of deviation Az2{t) at CO =0.4 and 0.6 a.u. is rather 
similar: in both cases electron 62 of atom B gets a momen- 
tum in the positive direction of the Z2-axis due to the Coulom- 
bic interaction with the laser-driven electron ei coming from 
atom A, which results eventually in the ionization of atom B 
in the positive direction of the Z2-axis. Note that a similar 
single-peak dynamics of electron 62 deviation Az2{t) occurs 
at all other laser carrier frequencies studied in this work (i.e. 
0.1 a.u. < CO <1.3 a.u.), except for the resonant frequency 
CO =0.5 a.u., when the Coulombic interaction of the two elec- 
trons is comparatively weak albeit durable (see Fig. |3] solid 
line). 

The time delay of the maxima for CO =0.4 and 0.6 a.u. ob- 
served in Fig. |3]finds a simple explanation in classical terms. 
With a kinetic energy of 0.1 a.u., which is available to the 
electron ei after excitation with co = 0.6 a.u. it can reach atom 



B within ^5.5 fs. Excitation with co = 0.4 a.u. requires a 
two photon transition for ionization and therefore the initial 
kinetic energy will be about 0.3 a.u. so that it takes only ^3 fs 
for electron 6\ to reach atom B. The delay between the max- 
ima of Az2(0 ill Fig- E] nicely agrees with this simple esti- 
mate. The resonant case CO = 0.5 a.u. is special since the 
created wave packet will be of Rydberg character. Such Ry- 
dberg wave packets have recently been studied for H atoms 
in ref. 18. There it was found that the maximum of the wave 
packet reaches '^150 a.u. on a time scale as long as 24 fs. 
Of course, in the present case this wave packet is subject to 
the Coulombic interaction with atom B and, therefore, it is 
strongly disturbed as compared with the single H-atom case. 

The resonant excitation of atom A at the laser carrier fre- 
quency CO =0.5 a.u. provides an example of an intuitiv6 case, 
when the ionization probability Ia of the laser-excited atom A 
is large compared to the ionization probability Ib of atom B 
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Fig. 5 Excitation of the extended H-H system by a laser pulse with the narrow spatial envelope of eqn (|5) acting only on electron (coordinate 
Zi) of atom A: electron probabilities P{zi) and P{z2) defined by eqn l llOt ) at two time moments, (a) - t =3 fs, CO =0.5 a.u.; (b) - t =10 fs, 
CO =0.5 a.u.; (c) - 1 =3 fs, CO =0.4 a.u.; (d) - / =10 fs, CO =0.4 a.u.. The other laser pulse parameters are: To =0.02 a.u. and tp =5 fs. Starting at 
±91 a.u. the wave function is absorbed by an imaginary optical potential. 



which is not excited by the laser field. The laser-driven dy- 
namics of the extended H-H system in the case of CO =0.5 a.u. 
is presented in the left panel of Fig. Hon a 10 fs timescale. 
For the sake of comparison, the right panel of Fig. 21 shows 
an example of a counter-intuitive case, CO =0.4 a.u., when the 
ionization probability Ib of atom B, which is not excited by 
the laser field, is substantially larger than the ionization prob- 
ability Ia of the laser excited atom A. 

From Figs. Ua) and (d) one can see that the 'atomic' en- 
ergies £'a(0 are controlled by the applied laser pulses: they 
increase in the first half of the pulses and decreases to the end 
of the pulses. In contrast, the 'atomic' energies E-Q{t) do not 
follow the applied laser pulses: they slowly increase in the sec- 
ond half of the pulses and after the end of the pulses. A similar 
behavior is found for 'atomic' ionizations Ia and Ib presented 
in Figs. Ub) and (e): while the laser-induced ionizations Ia 
rise fast in the second half of the pulses, ionization probabil- 
ities Ib start to rise only at the end of the pulses and sharply 
rise after the end of the pulses. In the counter-intuitive case 
presented in Fig. SJe), the ionization probability Ib is almost 
two times larger than Ia- 



Taking into account that the laser pulses with narrow spa- 
tial envelopes [Fig.[Tfb)] excite only electron ei of atom A in 
the vicinity of za = —50 a.u. and do not affect electron e2 
of atom B in the vicinity of zb — 50 a.u., we can assume that 
the energy transfer from A to B [Figs. 2Ia) and (d)] and, espe- 
cially, the time-delayed ionization of atom B in the positive di- 
rection of the Z2-axis [Figs.|4jb) and (e)], occur entirely due to 
electron-proton attraction Vep(zi,^) and the electron-electron 
repulsion Vee(zi,Z2) in the vicinity of atomB atZB =50 a.u.. 

In order to support these assumptions, we have calcu- 
lated the electronic probabilities P{z\) and P{z2), defined by 
eqn ( fTOb . during the laser pulse (f =3 fs) and after 5 fs free evo- 
lution of the H-H system after the end of the pulse (f =10 fs). 
The results obtained are shown in Fig. |5]for the two excita- 
tion conditions of Fig. |4] First, let us focus on the resonant 
excitation case shown in panels (a) and (b) of Fig. |5] Be- 
cause the laser pulse with the narrow spatial envelope excites 
only electron ey in the vicinity of za ~ —50 a.u., the probabil- 
ity distribution P{z\) during the pulse [at f = 3 fs. Fig. [Sja)] 
is broadened in comparison with the initial one [Fig. [Ha)] in 
both zi < ZA and zi > za directions. The laser-induced ex- 
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Fig. 6 Quantum dynamics of H-H excited from tlie entangled molecular singlet initial state by the narrowly shaped 5 fs laser pulse acting only 
on atom A in the vicinity of za = ^50 a.u.: the intuitive atomic ionization (left panel, to =0.5 a.u.) and the counter-intuitive atomic ionization 
(left panel, CO =0.4 a.u.). (a) and (c) - 'atomic' energies and defined by eqn J14b : (b) and (d) - total 'atomic' ionizations defined by 
eqns. l |15l l. 



tension of P{zi) into the domain of zi < za gives rise to the 
ionization Ia in the negative direction of the zi-axis. In con- 
trast, at zi > ZA, the laser-driven electron ei reaches the do- 
main zi > where it is attracted and accelerated by proton 
Pb, localized at zb = 50 a.u. [see Fig.|3b) for t = 10 fs]. 

The same qualitative behavior holds true for excitation with 
CO = 0.4 a.u. shown in Figs. |5jc) and (d). However, the details 
of the evolving wave packet are rather different giving rise to 
the different ionization probabilities. In particular we notice 
that due to the electron-proton Coulombic attraction, the prob- 
ability P{z\) has a local maximum at zi = Zb [Figs. |3c) and 
(d)], corresponding to laser induced electron transfer (LIFT) 
of ref. [ill. Simultaneously, when electron ei approaches the 
domain of zb = 50 a.u. and essentially localizes where the 
initial probability P{z2) of electron 62 has the global maxi- 
mum [Fig. m^)], the electron-electron Coulombic repulsion 
becomes very strong. Therefore, the probability P{z2) is con- 
siderably extended into the domain of Z2 > Zb, giving rise to 
ionization Ib in the positive direction of the Z2-axis. This effect 
is not so pronounced for resonant excitation where the wave 
packet describing electron ei is more delocaUzed in the region 
of atom B. 

Our numerical simulations performed with other laser car- 
rier frequencies have shown that the laser-driven dynamics of 
the electron probabilities P{zi) and P{z2) and their evolution 
after the end of the laser pulse are similar to those presented in 
Fig.|5]for CO =0.4 a.u.. We can conclude therefore, that the un- 



derlying physical mechanisms of the energy transfer from the 
laser-excited atom A to atom B are the Coulombic attraction 
of the laser-driven electron e i of atom A by the proton pB of 
atom B and a short-range Coulombic repulsion of electrons ei 
and 62 when their wave functions significantly overlap in the 
domain of atom B. 

3.2 Entangled Initial States 

In the previous sections we have have presented results ob- 
tained for the laser-driven dynamics of the extended H-H sys- 
tem of an "atomic origin" representing two distant H atoms. 
The initial state of such a system is an unentangled direct- 
product state. It is instructive to compare this case to two other 
alternatives: excitation of the H-H system from entangled ini- 
tial states, singlet and triplet, which represent an extended H-H 
system of a "molecular origin". 

For the sake of comparison to the intuitive and counter- 
intuitive cases of 'atomic' ionizations, presented in Fig.|4]for 
the unentangled atomic initial state, we performed numeri- 
cal simulations for the entangled molecular initial states, eq. 
( fT3] ), by making use of the same laser fields as for the unen- 
tangled atomic initial state: Tq =0.02 a.u., CO =0.5 a.u. and 
CO =0.4 a.u., with the narrow spatial envelope of the 5 fs laser 
pulse acting only on atom A. 

The laser-driven dynamics of extended H-H excited from 
the entangled molecular singlet state is presented in Fig. |6] 
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Notice that we have found that the respective resuhs for the 
initial triplet state are almost identical to those shown in 
Fig. |6]and therefore are not presented. We can see that the 
time-dependent 'atomic' energies £'a(0 and £'b(0 shown in 
Figs. |6fa) and |6jc) are quite different from those shown in 
Figs. Ua) and |4|d) for the unentangled atomic initial state. 
In contrast, the total time-dependent ionization probabilities 
7^°'"' and /g°'^', shown in Figs.|6jb) and|6jd), are very similar 
(yet not identical) to the time-dependent ionization probabili- 
ties /a and /b presented in Figs. Ub) andUe) for the case of 
the unentangled atomic initial state. In particular, the intuitive 
ionization (/a > /b) takes place in both cases at (O =0.5 a.u., 
while the counter-intuitive ionization occurs in both cases at 
CO =0.4 a.u.. 

From the comparison of the results presented in Figs. 2]and 
|6]the following observations can be made. 

(i) The overall energies transferred on a time scale of f = 
10 fs from atom A to atom B in the entangled molecular singlet 
and triplet states are almost identical to each other and similar 
to that transferred from A to B in the unentangled atomic state. 

(ii) In the case of the unentangled initial state, 'atomic' en- 
ergy starts to rise in the second half of the laser pulse and 
continues to rise after the end of the pulse, i.e., only when 
the laser-driven electron of atom A approaches the domain of 
atom B. In contrast, in the case of the entangled molecular ini- 
tial states, the 'atomic' energy of atom B, which is separated 
by 100 a.u. from the laser-excited atom A, is controlled by the 
laser pulse acting only on atom A: the energy (?) increases 
during the laser pulse and oscillates out-of -phase with Ep^{t) 
during the laser pulse and after the end of the pulse. Such 
an entangled behavior of distant atoms during coherent exci- 
tation of one of them is important evidence for the possibility 
of long-range quantum communication. 

(iii) The time-dependent ionization probabilities are almost 
identical in all three cases of the initial states of the extended 
H-H system. We can conclude therefore that the entangle- 
ment of the initial state of the extended H-H system does not 
change its ionization probability as compared to the unentan- 
gled atomic initial state. In particular, the intuitive atomic ion- 
ization remains intuitive and the counter-intuitive atomic ion- 
ization remains counter-intuitive. 

Note that similar observations have been made in our recent 
work— for the extended H-H system excited with a very differ- 
ent laser frequency of co = 1 .0 a.u.. We can conclude therefore 
that observations (i)-(iii) are quite generic. 

An interesting point concerns the question whether there is 
any difference at all in the laser-driven dynamics of the ex- 
tended H-H system excited from an entangled molecular ini- 
tial states, singlet and triplet? In order to give a quantitative 
answer to this question, we calculated the time-dependent in- 
tegrated probability Pee{t) to find two electrons at the same 
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Fig. 7 The time-dependent probabilities to find two electrons of the 
extended H-H system at the same point for the singlet and triplet 
initial states. The probabilities are defined by eqn ( I17t . the laser 
carrier frequency is to = 0.4 a.u.. 

point, 

p,,{t) = j duj dz\'y{R,z,z,t)\^. (17) 

The results obtained at the laser carrier frequency co = 0.4 a.u. 
are presented in Fig.|7](in the case of CO = 0.5 a.u., both prob- 
abilities are about four times smaller as compared to those 
shown in Fig. |7]i. First we notice that even for the singlet 
state the probabilities Pee{t) of eqn ( fTTI l are very small. The 
respective probabilities for the triplet initial state are always 
by about two orders of magnitude smaller in comparison to 
that of the singlet state. Therefore, the phase space available 
for electron-electron collisions in the triplet state is signifi- 
cantly smaller due to the antisymmetry of the electron wave 
function in the triplet state in comparison to that in the singlet 
state. Strictly speaking, Pee{t) should vanish for the initial 
triplet state and the finite but very small values shown in Fig. 
|2] reflect the numerical approximation. In order to study its 
influence we have performed a numerical experiment, where 
we set ^(/?,z,z,f) = by hand at every time step. The result 
for the 'atomic' energies and ionization dynamics are virtually 
indistinguishable form those shown in Fig. |6] 

The results presented in Fig. |7]provide also a visualization 
of the electron-electron collisions for the singlet initial state of 
the extended H-H system. Indeed, we can see that the proba- 
bility to find both electrons of H-H at the same point in space 
reveals strong pulsations in the second half of the 5 fs laser 
pulse and after the end of the pulse, i.e., when the laser-driven 
electrons of atom A approach the domain of atom B where 
they interact with electrons of atom B. What we see, in fact, is 
how the kinetic energy of the laser-driven electrons of atom A 
competes with the potential energy of thek repulsion by the 
electrons of atom B. 
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4 Conclusion 

In the present work we have studied numerically the non- 
Bom-Oppenheimer quantum dynamics of two distant H atoms 
(A and B) with an arbitrary but large initial intemuclear sepa- 
ration of 100 a.u. (5.29 nm), referred to as the extended H-H 
system. Both unentangled atomic states and entangled molec- 
ular states, singlet and triplet, were used as initial states in our 
simulations. 

First, the extended H-H system was assumed to be initially 
in an unentangled state, and only atom A was excited by the 
laser pulses with narrow spatial envelopes and various car- 
rier frequencies. Here, we have found that an efficient en- 
ergy transfer from atom A to atom B and ionization of atom B 
take place, being induced by a short-range electron-electron 
repulsion in the vicinity of atom B due to the preceding long- 
range transfer of the laser-driven electron ei from atom A en- 
hanced by the Coulombic attraction and acceleration of the 
laser-driven electron ei by proton of atom B. Furthermore, 
we observed that in the low-frequency domain CO < 0.6 a.u., 
EER-induced ionization of atom B depends very strongly on 
the carrier frequency of the laser pulse acting on atom A and 
can be both less efficient (an intuitive case) and more efficient 
(a counter-intuitive case) than the laser-induced ionization of 
atom A. 

Both intuitive (CO = 0.5 a.u.) and counter-intuitive (CO ~ 
0.4 a.u.) cases have been also studied for the entangled 
molecular states, singlet and triplet. We found that the time- 
dependent 'atomic' energies EA{t) and EB{t) obtained for the 
entangled molecular states are quite different from those ob- 
tained for the unentangled atomic initial state. In contrast, the 
total time-dependent ionization probabilities and /g"'"' 
for the entangled molecular initial states are very similar (yet 
not identical) to the time-dependent ionization probabilities 
I A and /b for the case of the unentangled atomic initial state. 
Moreover, both 'atomic' energies and ionization probabilities 
in the extended H-H system excited from the initial triplet state 
are almost identical to those of the singlet initial state. We 
found that the physical reason for the indistinguishability of 
the laser-driven dynamics of the extended H-H system excited 
from the initial singlet and triplet state is a very small prob- 
ability to find two electron of the extended H-H at the same 
point in space. 

We have shown in detail that in the case of the unentan- 
gled atomic initial state and a narrow spatial envelope of the 
applied laser field, when only the electron initially belonging 
to atom A is excited by the laser field, the physical mecha- 
nisms of the energy transfer from atom A to atom B and the 
ionization of atom B are as follows: (i) the Coulombic attrac- 
tion of the laser-driven electron ei of atom A by proton of 
atom B, resulting, in the case of CO 7^ 0.5 a.u., in the formation 
of a narrow local maximum of the electronic wave function 



in the vicinity of proton pB, and (ii) the short-range Coulomb 
repulsion of electrons ei and e2 in the vicinity of proton ps 
of atom B, where their wave functions strongly overlap. Due 
to the indistinguishability of the atomic ionization probabil- 
ities in extended H-H excited from the unentangled atomic 
state and from the entangled molecular states, both singlet and 
triplet, and on the basis of the results obtained in our recent 
work^, where the laser-driven dynamics of extended H-H has 
been studied for the unentangled atomic state and for the en- 
tangled molecular singlet state at CO = 1.0 a.u., we can con- 
clude that the physical mechanisms of the energy transfer from 
atom A to atom B and the ionization of atom B in the entan- 
gled states are basically the same as those in the unentangled 
atomic state. 

Note finally that the narrow spatial envelope of the laser 
pulse, exciting only one atom A of the extended H-H, used 
in the present work, is by no means unique, or obligatory for 
the long-range energy transfer to occur It has been shown in 
our recent work^ that very similar results are obtained when 
atom A is excited predeominantly at the edge of the Gaussian 
spatial envelope of the laser pulse. The latter case is more 
likely to arise in an experiment. 
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